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We present an interferometric and single dish study of small organic species 
toward Comets C/1995 01 (Hale-Bopp) and C/2002 T7 (LINEAR) using the 
BIMA interferometer at 3 mm and the ARO 12m telescope at 2 mm. For Comet 
Hale-Bopp, both the single-dish and interferometer observations of CH 3 0H indi- 
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cate an excitation temperature of 105±5 K and an average production rate ratio 
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Q(CH 3 OH)/Q(H 2 0)~1.3% at ~1 AU. Additionally, the aperture synthesis ob- 
servations of CH3OH suggest a distribution well described by a spherical outflow 
and no evidence of significant extended emission. Single-dish observations of 
CH 3 CN in Comet Hale-Bopp indicate an excitation temperature of 200±10 K 
and a production rate ratio of Q{CH 3 CN)/Q{H 2 O)~0. 017% at ~1 AU. The 
non-detection of a previously claimed transition of cometary (CH^OH^ toward 
Comet Hale-Bopp with the 12m telescope indicates a compact distribution of 
emission, D<9" (<8500 km). For the single-dish observations of Comet T7 LIN- 
EAR, we find an excitation temperature of CH3OH of 35±5 K and a CH3OH 
production rate ratio of Q(CH 3 OH)/Q(H 2 0)~1.5% at ~0.3 AU. Our data sup- 
port current chemical models that CH3OH, CH3CN and (CB^OH^ are parent 
nuclear species distributed into the coma via direct sublimation off cometary ices 
from the nucleus with no evidence of significant production in the outer coma. 

Subject headings: astrobiology - comets: individual (HALE BOPP (C/1995 01), 
LINEAR (C/2002 T7)) - molecular processes - techniques: interferometric - radio 
lines: solar system 



1. INTRODUCTION 

Over the last several years, complementary observations of comets with both single-dish 
telescopes and interferometric arrays, have shed new light on the abundance, production 
rate, distribution and formation of molecules in cometary comae. Recently, Milam et al. 
(2006) employed this method to investigate the origins of cometary formaldehyde (H2CO) 
by comparing observations with both the Arizona Radio Observatory (ARO) 12m 12 l single- 
dish telescope and the Berkeley-Illinois-Maryland Association (BIMA) array 13 !. From this 
investigation toward Comet C/1995 01 (Hale-Bopp), the authors clearly showed that the 
distribution and abundance of H 2 CO must be coming from a source other than the comet nu- 
cleus. The extended source of H 2 CO may be due to an enhanced production of this molecule, 
(~4x larger) from grains consisting of silicates and organic material (Milam et al. 2006), or 
possibly from the photodissociation of a larger parent species such as methanol (CH 3 OH) 
(Hudson, Moore & Cook 2005). These complementary observations, however, are very dif- 



12 The Kitt Peak 12m telescope is currently operated by Arizona Radio Observatory (ARO), Steward 
Observatory, University of Arizona. 

13 Operated by the University of California, Berkeley, the University of Illinois, and the University of 
Maryland with support from the National Science Foundation. 
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ficult to obtain, as most of the molecular observations toward comets are performed using 
solely single-dish telescopes or interferometric arrays, depending on the telescope available 
at the time as the comet approaches the inner solar system and perigee. 

Multiple facility collaborations for cometary observations are necessary if we are to fur- 
ther understand the formation, excitation and distribution of organic cometary material. 
Recent cometary models of the abundance of organic material suggest Oort Cloud comets 
may have had an origin in the giant planet region, where they were thermally and chemically 
processed, before being ejected into the Oort Cloud (Weissman 1999). Furthermore, obser- 
vations comparing long- and short-period comets suggest distinct classes - organic "rich", 
"normal", and "depleted" (Bonev et al. 2008; Mumma et al. 2003). However, the current 
sample size is too limited to accurately characterize each class of comets. It should also 
be noted that many of the comets described as organic "rich" have undergone breakup or 
fragmentation events, thus exposing more internal, unprocessed material than an object just 
entering the inner solar system (Bonev et al. 2008). So, the question persists, how does the 
classification of organic "rich" , "normal" or "depleted" hold with regards to comets that are 
for example, dust rich versus dust poor, having undergone a recent breakup, or bear high 
abundances of hypervolatile species? 

The organic species examined in this study, methanol (CH3OH), methyl cyanide (CH3CN) 
and ethylene glycol ([CH 2 OH] 2 ) are believed to be parent species directly sublimated off the 
cometary ices in the nucleus. Since the first reported single-dish detections of cometary 
CH3OH towards Comet Brorsen- Met calf (1989 o) and Comet Austin (1989 cl) (Colom et al. 
1990), CH3OH remains the most abundant large (>5 atoms) organic (containing a C-atom) 
molecule measured in cometary ices. In a recent review by Charnley & Rodgers (2008), 
CH3OH has a measured abundance range relative to H 2 between 1-4%. However, aper- 
ture synthesis observations of cometary CH 3 OH are very rare and no high spatial resolution 
data of this molecule currently exists. Similarly, single-dish observations of CH 3 CN and 
(CH 2 OH) 2 also suggest a purely nuclear origin (Biver et al. 2002, Crovisier et al. 2004). 
CH3CN, as well as HC3N, may be parent species that would increase the abundance of the 
CN radical as well as C 2 in the coma. Yet, current observations do not preclude the possibil- 
ity of CH 3 CN having its own parent species (Bockelee-Morvan et al. 1987). In addition, the 
high abundance of (CH 2 OH) 2 observed in Comet Hale-Bopp at high spatial resolution also 
strongly suggests a nuclear origin. To date, complementary observations of the same transi- 
tions of either CH3CN or (CH 2 OH) 2 with another facility were not available to follow-up on 
the claimed distribution of these species in the coma and their possible formation pathways. 

In this paper, we report on complementary observations of cometary CH3OH, CH3CN 
and (CH 2 OH) 2 taken with the BIMA array near Hat Creek, CA and the ARO 12m telescope 
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on Kitt Peak, AZ toward Comets C/1995 01 (Hale-Bopp) and C/2002 T7 (LINEAR). From 
these observations, we continued to investigate the abundance, production rate, distribution 
and formation of molecules in comets by directly comparing the measurements of these 
organic species with different telescopes. Our data supports the hypothesis that CH 3 OH, 
CH3CN and (CH20H)2 are primarily parent species, generated by direct sublimation of 
cometary ices contained in the nucleus. 



2. OBSERVATIONS 

2.1. BIMA Array 

Our observations of Comet Hale-Bopp with the BIMA array constitute one of the earli- 
est detections of cometary CH 3 OH with an interferometer. These measurements were taken 
using the "soft" C configuration of the BIMA array from 1997 March 24 to April 03. Data 
were acquired in the interferometric (cross-correlation) mode with nine antennas. The min- 
imum baseline for these observations was ~15 m and the maximum baseline was ~139 m. 
The average resulting full width half power synthesized beam size was ~8" x 10" for the ob- 
servations. The quasar 0102+584 was used to calibrate the antenna based gains and the 
absolute amplitude calibration of this source was based on planetary observations and is ac- 
curate to within ~20%. The spectral windows containing the transitions had a bandwidth of 
25 MHz and were divided into 256 channels for a spectral resolution of 0.1 MHz per channel. 
However, to increase the signal-to-noise ratio in each window, the data were averaged over 
two channels giving an effective spectral resolution of 0.2 MHz per channel. All data were 
first corrected to JPL ephemeris reference orbit 139. The data were then combined and im- 
aged using the MIRIAD software package (Sault, Teuben & Wright 1995). To include all the 
data from the comet with multiple tracks, the data were inverted in u-v space. The resulting 
synthesized beamsize was determined from the combined dataset. The uvdata from multiple 
tracks were combined to make images. Table 1 summarizes the observational parameters of 
the BIMA array observations. 



2.2. 12m Single Dish 

Observations at 2mm of CH 3 OH, CH 3 CN and (CH 2 OH) 2 toward Comet Hale-Bopp 
were taken during three observing runs on 1997 Mar 10, 11, and 20 using the ARO (at 
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the time, NRAQ 14 !) 12m telescope. Observations of CH3OH toward Comet LINEAR were 
conducted on 2004 May 20 also using the ARO 12m telescope. Dual-channel SIS mixers 
were used for the 2 mm observations, operated in single- sideband mode with ~20 dB im- 
age rejection. The backends employed for the presented observations were individual filter 
banks with 500 kHz resolution and the millimeter autocorrelator (MAC) with a resolution 
of 781 kHz. The spectral temperature scale was determined by the chopper-wheel method, 
corrected for forward spillover losses, given in terms of Tp (K). The radiation temperature, 
Tr, is then derived from the corrected beam efficiency, rj c , where Tr = T^/t] c . A two body 
ephemeris program was used to determine the comet's position using the orbital elements 
provided by D. Yeomans (1996, private communication) of JPL for Comet Hale-Bopp and 
JPL ephemeris reference orbit 69 was used for the Comet LINEAR observations. Focus and 
positional accuracy were checked periodically on nearby planets or maser sources. Data were 
taken in the position-switched mode with the off position 30' west in azimuth. Observing 
frequencies, dates, comet geocentric and heliocentric distances, beam size, the diameter of 
the projected beam size on the comet, and the corrected main beam efficiency (rj c ) at the 
times of observations are listed in Table 2. 

Table 3 lists the transitions and molecular line parameters of CH 3 OH, CH 3 CN and 
(CH20H)2- All spectral line data were taken from the Spectral Line Atlas of Interstellar 
Molecules (SLAIM) available at www.splatalogue.net (Remijan et al. 2008) and the Cologne 
Database for Molecular Spectroscopy (Miiller et al. 2005). 



3. RESULTS 

3.1. BIMA Array 

Figure 1 shows the images and spectra of CH 3 OH taken toward Comet Hale-Bopp with 
the BIMA array. Figure 1(a) shows the image of the J(K a ,K c )=2(0,2)-l(0,l) E transition 
of CH3OH averaged over 2 days, 1997 March 27 and 31. Notice that the CH3OH emission 
is slightly larger than the synthesized beam of the BIMA array (located at the bottom 
left of each image) and follows the general position angle of the synthesized beam. The 
projected direction to the sun is shown by the line segment. The coordinates are given in 
offset arcseconds centered on the comet nucleus. Figure 1(b) shows the spectra of 4 CH3OH 
transitions near 96.7 GHz (Table 3) including the J(K a ,K c )=2(0,2)-l(0,l) E transition. The 



14 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by the Associated Universities, Inc. 
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spectrum has been hanning smoothed over 3 channels. The rms noise level is ~0.1 Jy beam" 1 
(indicated by the vertical bar at the left of the spectrum). The spectral line labels correspond 
to the rest frequency located at the top left of the panel for a cometocentric velocity of 
km s _1 , denoted by a dashed line. 

Figure 1(c) shows the image of the J(K a ,K c )=3(l,3)-4(0,4) A+ transition of CH 3 OH 
averaged over 3 days, 1997 March 26 and April 2-3. Figure 1(d) shows the spectrum of this 
CH 3 OH transition near 107.0 GHz (Table 3). Note that the two CH 3 OH emission contours 
shown in figures (a) and (c) are centered on the predicted location of the comet nucleus 
with no apparent offset. The images are similar to each other and similar to the average of 
the March HCN BIMA array contours given by Veal et al. (2000) and the CO BIMA array 
contours given by Milam et al. (2006). In addition, there is no pronounced evidence for the 
jets or significant extended enhancement of CH3OH beyond the 3<r detection level (§4.2). 
Furthermore, from figures 1(b) and (d), there is no apparent second velocity component of 
CH3OH nor is the emission offset from the cometocentric rest velocity of km s _1 . 



3.2. 12m Single Dish 

Figure 2 shows the spectra of the five A ~ 2mm transitions of CH 3 OH detected toward 
Comet Hale-Bopp with the 12m telescope (Table 3). The UT date and spectral line quantum 
numbers of each observed transition are located at the top of each panel. As in the BIMA 
spectra, the spectral line labels correspond to the rest frequency of that transition for a 
cometocentric velocity of km s" 1 (dashed line). The CH 3 OH detections reported by the 
12m data are all E-type transitions. Similarly, Figure 3 shows the spectrum of a cluster of 
E-type transitions of CH 3 OH centered around 157.3 GHz observed toward Comet LINEAR 
with the ARO 12m telescope. 

Finally, Figure 4 shows the spectrum taken at the frequency of the J = 8 — 7 (K=0 
through 7) transitions of the symmetric top species methyl cyanide (CH 3 CN) toward Comet 
Hale-Bopp using the 12m. The K=0 through 6 components are clearly present in these data. 
CH3CN and other symmetric tops have properties that make them ideal probes of physical 
conditions of astronomical environments (see e.g. Araya et al. 2005; Remijan et al. 2004; 
Hofner et al. 1996 ). As with CH 3 OH, there is no apparent additional velocity components 
of CH 3 CN or offset from the cometocentric rest velocity of km s _1 . 

Table 4 lists the molecular species that were searched for toward Comets Hale-Bopp and 
LINEAR at or near the cometocentric velocity (0 km s" 1 ). Least square Gaussian fits were 
made for each spectral line in order to obtain the peak intensities and line-widths for the 
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detected transitions. For those transitions not detected, the 1 a rms noise level is given. The 
spectroscopic parameters used in determining the column densities and production rates are 
found in Table 3 and the formalism for calculating column densities and production rates is 
described in 84.1. 



4. ANALYSIS AND DISCUSSION 

4.1. EXCITATION, ABUNDANCES AND PRODUCTION RATES 

Determining the best temperature, column density, and production rates of cometary 
molecules are important because they are essential to understanding the formation of cometary 
species. This includes the role of grains in molecule formation, the relative impact of warm 
versus cold gas phase chemistry, or the importance of photodissociation in the production of 
daughter species in the coma. In order to determine the temperatures and column densities 
of the region of the coma that contains the CH3OH and CH3CN emission, we assume that 
each region has uniform physical conditions, that the populations of the energy levels can 
be characterized by a Boltzmann distribution and finally, that the emission is optically thin. 
The rotational temperature diagram method was not employed for this analysis because even 
small errors in the intensities and widths of inherently weak spectral lines can cause large 
errors in determining the temperature and total column density (see e.g. Snyder et al. 2005). 
The complete formalism for these calculations using the BIMA array data can be found in 
Remijan et al. (2006) and for the 12m data, in Milam et al. (2006). 

Figure 1 shows the results of the model fit (blue trace) applied to the CH3OH BIMA 
array data. The array data are well fit by an excitation temperature of T ex =105±5 K and 
a total beam-averaged column density of N^=9 ± 1 x 10 14 cm~ 2 assuming that the CH 3 OH 
emission is extended beyond the 10."33x8."46 average synthesized beam of the BIMA array 
(i.e. the beam filling factor, f=l). Similarly, Figure 2 shows the results of the model fit 
applied to the CH3OH ARO 12m data with the same excitation temperature and a slightly 
lower column density of N^=2.4 ± 0.3 x 10 14 cm -2 . The average total column density from 
these two facilities is N^=5.7±0.7 x 10 14 cm -2 . As the figure shows, the modeling reproduces 
the observations. Assuming CH 3 OH is a parent species, production rates were determined 
from a Monte Carlo model for a nuclear source (see Combi & Smyth, 1988). This model 
traces the trajectories of molecules, within the telescope beam, ejected from the comet 
surface. The observed column density is then matched for an output molecular production 
rate, Q (s _1 ). Input parameters include the lifetime of the molecule, from Huebner et al. 
(1992), scaled by R^, the outflow velocity, telescope beam size, and the observed column 
density. Assuming a nuclear origin for CH3OH and using this Monte Carlo model, we find 



an average CH3OH production rate from both instruments of Q(C H^OH)^l. 3 xlO 29 s l or 
Q(CH 3 OH)/Q(H 2 0)~1.3% at -1 AU. This production rate ratio assumes a Q(H 2 O) = 10 31 
s _1 (Crovisier et al. 2004). 

In Figure 2, there is a slight discrepancy between the observed intensities and model 
predictions for the CH3OH lines at 151.86 GHz. The model predicts a higher line intensity 
than observed and the model underestimates the line intensity at 148.11 GHz and 145.76 
GHz. However, given the S/N ratio of all the passbands where these lines are observed, 
the model does a very good job in predicting the expected line intensity of each CH3OH 
transition observed with the 12m telescope. The same can be said for the array data. The 
only discrepancy is the 96.7445 GHz, J(K a ,K c )=2(0,2)-l(0,l) E CH 3 OH transition (Figure 
1); the intensity is somewhat underestimated by our model fit. 

Figure 4 shows the best model fit to the J(K) = 8(K) — 7(K) (K=0 through 6) 
transitions of methyl cyanide (CH3CN) toward Comet Hale-Bopp. Once again, the blue 
trace shows the model predictions for the CH3CN data which is best fit by an excitation 
temperature, T e;c =200±10 K and a total beam-averaged column density, N^=2.6 ± 0.3 x 
10 12 cm -2 . In this case, the predicted K = 7 line intensity is below the 1 o detection 
limit, which indicates that the emission feature near the K = 7 rest frequency may not 
be real. Assuming a nuclear origin for CH3CN and using a Monte Carlo model, we find a 
production rate of Q(C# 3 CA0~1.7xl0 27 s" 1 or Q(CH 3 CN)/Q(H 2 O)~0.017% at ~1 AU. 
This production rate is similar to what Biver et al. (1997) found, ~10 27 s _1 at 1 AU, in their 
study. 

These 12 m data overlap in frequency with the (CH20H)2 IRAM 30m study of Crovisier 
et al. (2004). In the IRAM dataset, the claimed detection of the J(K a ,K c ,u)=15(7,9,0)- 
14(7,8,1) and 15(7,8,0)-14(7,7,1) transitions at 147.13 GHz are nearly coincident in frequency 
and show up as one feature given the spectral resolution. Moreover, the (CH 2 OH) 2 transitions 
have an overall line intensity equal to that of the J(K) = 8(4) — 7(4) transition of CH3CN 
detected with the IRAM 30m. From the ARO 12m data, we identify the same spectral 
features of CH 3 CN but do not detect any emission from the two transitions of (CH 2 OH) 2 
below the la limit. From the CH 3 CN, CH 3 OH and (CH 2 OH) 2 12m data, it is clear we are 
probing different physical environments in the coma of Comet Hale-Bopp. A discussion of 
the distribution of these molecular species is in §4.2. 

Figure 3 shows the best model fit (blue trace) to the cluster of E-type transitions of 
CH3OH centered around 157.26 GHz toward Comet LINEAR with the ARO 12m tele- 
scope. There have been very few molecular line observations toward Comet LINEAR. 
The most complete study of molecular species toward this object was by Remijan et al. 
(2006) using the BIMA array. These authors detected the J(K a ,K c )=3(l,3)-4(0,4) A+ tran- 
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sition of CH3OH with a synthesized beam size of 22. "3 x 16. "9. Using the proposed ex- 
citation temperature of 115 K (DiSanti et al. 2004; Magee-Sauer et al. 2004; Kiippers et 
al. 2004), Remijan et al. determined a total beam-averaged column density of CH 3 OH of 
Ny=1.4 ± 0.3 x 10 14 cm" 2 . However, as shown in figure 3, the 12m CH 3 OH data are 
best fit by an excitation temperature of T e:r =35±5 K and a total beam-averaged col- 
umn density of Nr=2.2 ± 0.3 x 10 13 cm" 2 (blue trace), implying a production rate of 
Q(CH 3 OH)~2.0xl0 27 s" 1 or Q(CH 3 OH)/Q(H 2 0)~1.5% at 0.3 AU. This production rate 
ratio assumes a Q(H 2 O)=1.3xl0 29 s" 1 (Lecacheux et al. 2004). To reconcile the differing 
column density determined by Remijan et al. (2006), we used the best model fit parameters 
of temperature and column density from the 12m data to predict the line intensity of the 
J(K a ,K c )=3(l,3)-4(0,4) A + transition detected by the BIMA array. Figure 5 shows the result 
of this analysis (5(a) shows the hanning smoothed data and 5(b), the raw data as shown 
in figure 1 of Remijan et al. (2006)). The fit is excellent. From using the ARO data it 
is quite clear that the assumption of an excitation temperature of T ex =115 K for CH3OH 
was incorrect and the determined column density of CH3OH toward Comet LINEAR was 
therefore overestimated by almost an order of magnitude. 



4.2. DISTRIBUTION OF MOLECULAR SPECIES 

We present two complementary sets of data from Comet Hale-Bopp that give great in- 
sights into the distribution of CH 3 OH, CH 3 CN and (CH 2 OH) 2 . First, it is quite clear from 
the data presented in Figures 1 and 2 that despite the drastic change in the telescope beam 
sizes between the array and the 12m, we are still sampling a region of the coma that is ther- 
malized to ~105 K. The ratio of areas subtended by the 12m beam, compared to the BIMA 
array synthesized beam, is ~21:1 which corresponds to a change in the physical sampled size 
scale from D~40000 km to D~8500 km. A direct comparison of the column densities from 
the 12m and BIMA array suggest an enhancement of the CH3OH abundance at smaller coma 
radii. Furthermore, Column 5 in table 3 lists the Einstein A coefficients for each transition 
of CH3OH, CH 3 CN and (CH 2 OH) 2 . Assuming a collision cross section of ~10 -15 cm 2 and 
a temperature of 105 K, the calculated critical densities of the detected CH3OH transitions 
range from (0.7-3.2) xlO 5 cm" 3 . Our observations encompass a region with a H 2 density of 
^lo^ 5 - 7 " 7 - 5 ) cm" 3 (Lovell et al. 2004). Therefore, there is sufficient collisional excitation to 
populate the CH 3 OH levels, and our observations are sampling molecular abundances and 
not excitation effects. 

Figure 6(a) once again shows the distribution of the J(K a ,K c )=3(l,3)-4(0,4) A + transi- 
tion of CH3OH around Comet Hale-Bopp. Figure 6(b) is a model fit to the data shown in (a). 
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The model was generated assuming a spherical Haser model of parent species sublimating 
off the comet's nucleus. From the data in (a), the scale length of the distribution of CH 3 OH 
around the nucleus is measured to be r~10". Using this scale length and the MIRIAD task 
IMGEN, a Haser model distribution of CH 3 OH was generated which was then convolved 
with the synthesized beam of the BIMA array (shown at the bottom left of each panel in 
figure 6). The residuals obtained by subtracting the model from the data are shown in figure 
6(c). While it appears that there may be a more extended distribution of CH 3 OH detected 
by the BIMA array based on the residual data, it is important to note that the contour levels 
shown in (a) and (b) are 3, 4, 5 and 6 a (cr=0.05 Jy beam™ 1 ), whereas the data shown in 
(c) are 1 and 1.5 a. Although there is some indirect evidence for an extended distribution of 
CH3OH with D = 40000 km from the 12 m detections of CH3OH, the data presented in (c) 
indicate that any extended emission is is completely resolved out by the array. The array 
would be sensitive to extended emission of CH 3 OH in the outer coma only if the distribution 
was "clumpy", as is seen in H2CO (Milam et al. 2006). 

It appears that the CH3CN is tracing higher temperature and higher density gas than 
CH3OH. Because CH3CN has a larger dipole moment than CH3OH, its measured excitation 
should be lower, given identical conditions. However, if the overall distribution of CH 3 CN 
is closer to the nucleus than CH 3 OH, we would expect to measure a higher excitation tem- 
perature (it is well known that the temperature of the coma increases with decreasing coma 
radii (Biver et al. 1999)). This distribution would also imply that CH3CN is a parent species 
(§4.3). Given the difference in beam sizes between the 12m and 30m telescopes, we can get 
an estimate of the size of the emitting region based on the relative beam sizes using the beam 
filling factor in the expression for column density. Comparing the relative line intensities 
between the K=3 and 4 transitions of CH3CN from the IRAM 30m data and what is mea- 
sured from the 12m, the 30m emission features from CH 3 CN are ~5-10x stronger than what 
is observed from the 12m. This difference implies either 1) an enhancement of the CH 3 CN 
column density at smaller coma radii or 2) the excitation of the transitions we observed of 
CH3CN takes place at a slightly smaller radii where the density is higher and the measure- 
ments taken with the 12m are suffering from some beam dilution. Assuming the H2O density 
profile of Lovell et al. (2004), the maximum coma radius that lies below the CH 3 CN critical 
density range of (1.2-5.3) xlO 6 cm -3 is ~10000 km. This explains why the 30m CH 3 CN line 
strengths are higher than those from the 12m. Finally, from the non-detection of (CH20H)2 
in the 12m observations, it is also apparent that the distribution of (CH20H)2 must also 
be compact in the cometary coma. Assuming the detected transition line strength of ~0.09 
K in the 30m spectrum is from (CH 2 OH) 2 and a la line intensity of ~0.01 K in the 12m 
passband, the emission from (CH 2 OH) 2 must be on the order of D<9" corresponding to a 
physical size of D<8500 km. High sensitivity, high resolution interferometer observations 



-11- 



are necessary to confirm this distribution. 



4.3. FORMATION OF COMETARY METHANOL AND OTHER 

MOLECULAR SPECIES 

From the high resolution methanol observations of the BIMA array, it is clear that 
the CH3OH is either sublimating directly off the cometary ices contained in the nucleus 
or is formed very deep in cometary coma. If CH3OH originates in the ices, this molecules 
may be a remnant from the formation of the presolar nebula and hence the interstellar 
medium. Currently, there are two accepted formation pathways to the interstellar production 
of CH3OH, including the radiative association of CH^ + H 2 —> CH 3 OH^", followed by 
recombination with an electron to produce CH3OH and H (Herbst, 1985). Grain surface 
reactions that involve the repeated hydrogenation of CO — >■ H + CO — > HCO + H — > 
HCHO -> HCHOH are thought to eventually lead to CH 3 OH (Hiraoka et al. 1994). Our 
observations show there is no significant enhancement in the production of CH 3 OH in the 
outer coma as the comet enters the inner solar system. 

The same argument can be applied to CH3CN. The accepted route to the formation 
of CH3CN in interstellar environments is CH^ + HCN, a collision complex is formed that 
equilibrates to CH 3 CNH + + v (Herbst, 1985); then CH 3 CNH + combines with an electron to 
form CH 3 CN + H. Our observations of CH 3 CN show that it is present close to the cometary 
nucleus, indicating direct sublimation off cometary ices. However, the excitation of the 
transitions we observed of CH3CN takes place at a smaller radii where the density is higher, 
and this molecule may be present as well at larger radii. However, the true spatial scale and 
scale length of CH 3 CN needs to be verified by interferometric observations. 

Finally, the non-detection of (CH20H)2 with the ARO 12m gives insight into the for- 
mation of (CH 2 OH) 2 as well as the possibility of detecting CH 2 OHCHO. If the transition 
detected by the 30m is indeed from cometary (CH 2 OH) 2 , in §4.2, we found the predicted 
distribution of (CH 2 OH) 2 toward Comet Hale-Bopp was calculated to be D<9" (<8500 
km), again suggesting direct sublimation off cometary ices or an enhanced production in 
the inner coma. Presumably, as with the formation of CH3OH and CH3CN in interstellar 
environments, (CH 2 OH) 2 would be formed in the ISM and then seeded into cometary ices. 
Currently, there is no accepted gas phase interstellar formation mechanism that can lead 
to the production of (CH 2 OH) 2 . However, there are several formation pathways using sur- 
face chemistry which may occur on icy grain mantles. For example, to form (CH 2 OH) 2 and 
CH 2 OHCHO, formaldehyde (H 2 CO) produces CH 2 OHCHO in an aqueous Formose reaction, 
which can then be hydrogenated to form (CH 2 OH) 2 (e.g., see Walker 1975). Charnley (2001) 
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predicts that (CH20H)2 and CH2OHCHO could be formed by direct hydrogenation reac- 
tions, starting from ketene (CH 2 CO) on grain surfaces. Finally, Hudson & Moore (2000) 
showed that proton irradiated CH 3 OH on icy interstellar grain mantles can lead to the for- 
mation of (CH 2 OH) 2 . Thus, there appear to be several ways to form and embed (CH 2 OH) 2 
into cometary ices. 



5. CONCLUSIONS 

We presented an interferometric and single dish study of cometary molecules toward 
comets C/1995 01 (Hale-Bopp) and C/2002 T7 (LINEAR) using the Berkeley-Illinois- 
Maryland Association (BIMA) interferometer at 3 mm and the Arizona Radio Observatory 
(ARO) 12m telescope at 2 mm. The overall conclusions from our analysis of these data are: 

1. The CH3OH Hale-Bopp data are well fit by an excitation temperature of T e;c =105±5 
K and a total beam-averaged column density of Np=5.7 ± 0.7 x 10 14 cm~ 2 . Assum- 
ing a nuclear origin for CH3OH and using a Monte Carlo model, we find an aver- 
age CH3OH production rate from both instruments of Q (CH 3 OH)~l. 3 xlO 29 s _1 or 
Q(CH 3 OH)/Q(H 2 0)~1.3% at - 1 AU. 

2. The CH3OH ARO 12m LINEAR data centered around 157.26 GHz are best fit by 
an excitation temperature of T e;c =35±5 K and a total beam-averaged column density 
of N T =2.2±0.3xl0 13 cm" 2 and thus, a production rate Q(CH 3 OH)~2.0xl0 27 s" 1 or 
Q(CH 3 OH)/Q(H 2 0)=1.5% at 0.3 AU. 

3. From the combination of the single-dish and aperture synthesis observations of CH 3 OH, 
we find the distribution of CH 3 OH toward Comet Hale-Bopp is well described by a 
spherical outflow with an increase in column density closer to the cometary nucleus. 
The data presented in the array images show no evidence of significant enhanced pro- 
duction of CH3OH in the extended coma or from jets as any extended emission is 
completely resolved out by the array. 

4. The CH3CN 12m Hale-Bopp data is best fit by an excitation temperature, T ex =200±10 
K and a total beam-averaged column density, Nr=2.6±0.3xl0 12 cm~ 2 . Assuming a 
nuclear origin for CH3CN and using a Monte Carlo model, we find an average CH3CN 
production rate of Q(CH 3 CN)~1.7xl0 27 s" 1 or Q(CH 3 CN)/Q(H 2 O)~0 .017% at ~ 
1 AU. A comparison between the single-dish observations from the ARO 12m and 
the IRAM 30m of CH3CN suggest that the ARO observations are beam diluted. The 
excitation of the transitions observed of CH3CN takes place at a smaller radii where 
the density is higher and suggest a nuclear origin of CH 3 CN. 
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5. The non-detection of a previously claimed transition of cometary (CH20H)2 toward 
Comet Hale-Bopp with the ARO 12m telescope indicates a compact distribution of 
emission on the order of <9" (<8500 km). This supports the hypothesis that the 
cometary production of (CH 2 OH) 2 is direct sublimation off cometary ices from the 
nucleus. 
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Table 1. CH 3 OH Observational Parameters - BIMA array 



Frequency 


Observation 


A 


r 


beam 


chan rms 


(GHz) 


date 


(AU) 


(AU) 


("x") 


(Jy beam -1 ) 


96.741 


1997 Mar 27 


1.323 


0.919 


10.8x8.2 


0.21 




1997 Mar 31 


1.344 


0.914 


9.8x7.8 


0.20 


107.014 


1997 Mar 26 


1.320 


0.921 


19.3x6.8 


0.43 




1997 Apr 02 


1.359 


0.914 


11.7x6.6 


0.15 




1997 Apr 03 


1.367 


0.915 


16.0x6.4 


0.30 
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Table 2. Observational Parameters - ARO 12m 



Frequency 


Observation 


A 


r 


beam 


D 


Vc 


molecule 


(GHz) 


date 


(AU) 


(AU) 


(") 


(km) 






Comet Hale-Bopp 


145.766 


1997 Mar 20.94 


1.318 


0.940 


43 


41105 


0.76 


CH 3 OH 


147.105 a 


1997 Mar 20.66 


1.318 


0.940 


43 


41105 


0.76 


CH 3 CN 


148.112 


1997 Mar 11.89 


1.368 


0.989 


42 


41672 


0.76 


CH3OH 


150.142 


1997 Mar 10.71 


1.377 


0.996 


42 


41946 


0.75 


CH3OH 


151.860 


1997 Mar 20.92 


1.318 


0.940 


41 


39193 


0.74 


CH3OH 


157.049 


1997 Mar 21.54 


1.316 


0.936 


40 


38179 


0.72 


CH3OH 




Comet LINEAR 










157.261 b 


2004 May 21.80 


0.284 


0.865 


40 


8239 


0.72 


CH3OH 



a Average of observing frequencies between the K—l and K=0 lines of 
CH 3 CN. 

b Average of observing frequencies of CH 3 OH lines detected in the same spec- 
tral passband. 
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Table 3. Molecular Line Parameters 



Molecule 


Transition 


Frequency 


< Sij/j, 2 > 


A*,3 


E„, 






(MHz) 


(Debye 2 ) 


(s- 1 ) 


(K) 


CH 3 OH a 


J(K a ,K c )=2(-l,2)-l(-l,l) E 


96,739.363(3) 


1.2 


2.55xl0~ 6 


12.6 




2(0,2)-l(0,l) A+ 


96,741.377(3) 


1.6 


3.41xl0~ 6 


7.0 




2(0,2)-l(0,l) E 


96,744.549(3) 


1.6 


3.41xl0~ 6 


20.1 




2(1,1)-1(1,0) E 


96,755.507(3) 


1.2 


2.62xl0~ 6 


28.0 




3(l,3)-4(0,4) A+ 


107,013.770(13) 


3.0 


6.13xl0~ 6 


28.4 




16(0,16)-16(-1,16) E 


145,766.163(27) 


7.0 


7.60xl0~ 6 


327.9 




15(0,15)-15(-1,15) E 


148,111.919(24) 


7.3 


8.80xl0~ 6 


290.9 




14(0,14)-14(-1,14) E 


150,141.593(22) 


7.5 


10.2xl0~ 6 


256.3 




13(0,13)-13(-1,13) E 


151,860.170(20) 


7.6 


11.5xl0~ 6 


223.9 




6(0,6)-6(-l,6) E 


157,048.586(13) 


5.7 


1.96xl0~ 6 


61.9 




4(0,4)-4(-l,4) E 


157,246.041(14) 


4.2 


21.0xl0~ 6 


36.4 




1(0,1)-1(-1,1) E 


157,270.818(15) 


1.5 


22.1xl0~ 6 


15.5 




3(0,3)-3(-l,3) E 


157,272.320(14) 


3.3 


21.5xl0~ 6 


27.1 




2(0,2)-2(-l,2) E 


157,276.004(14) 


2.4 


21.8xl0~ 6 


20.1 


CH 3 CN a 


J(K)=8(6)-7(6) 


147,072.608(2) 


53.2 


0.58xl0~ 4 


289.0 




8(5)-7(5) 


147,103.741(1) 


74.2 


1.62xl0~ 4 


210.5 




8(4)-7(4) 


147,129.232(1) 


91.3 


1.99xl0~ 4 


146.2 




8(3)-7(3) 


147,149.068(1) 


104.6 


1.14X10" 4 


96.1 




8(2)-7(2) 


147,163.243(1) 


114.1 


2.49xl0~ 4 


60.4 




8(1)-7(1) 


147,171.750(1) 


119.8 


2.61xl0~ 4 


38.9 




8(0)-7(0) 


147,174.587(1) 


121.7 


2.66xl0' 4 


31.8 


(CH 2 OH) 2 b 


J(K a ,K c ,u)=15(7,9,0)-14(7,8,l) 


147,131.814(36) 


2.5 


5.99xl0~ 5 


83.1 




15(7,8,0)-14(7,7,1) 


147,132.412(36) 


2.7 


5.99xl0~ 5 


83.1 



a Molecular line parameters of CH3OH and CH3CN taken from SLAIM at www.splatalogue.net. 
Q r =1.2Tf for CH3OH and Q r =3.9Tf for CH 3 CN. 

b Molecular line parameters of (CHjOH^ taken from CDMS (Miiller et al. 2005). The energy 



levels are noted as J(K a ,K c 
(Christen & Miiller 2003). 



where v is the quantum number associated with OH tunneling 
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Table 4. Comet Hale-Bopp and Comet LINEAR Molecular Line Identifications 

Species Telescope v Transition ^line & ^ v ^T Qp Qx/Qi^oO 

(MHz) (kms -1 ) (cm -2 ) (s _1 ) 

Comet Hale-Bopp 

CH 3 OH BIMA 96,739.4 2(-l,2)-l(-l,l) E 0.31(5) 2.4(2) 9±lxl0 14 l.l±0.3xl0 2a ~l.lxlO _21 

96.741.4 2(0,2)-l(0,l) A+ 0.52(5) 2.0(2) 

96.744.5 2(0,2)-1(0,1) E 0.60(5) 2.1(2) 
96,755.5 2(1,1)-1(1,0) E 0.29(5) 3.2(2) 

107.013.8 3(l,3)-4(0,4) A+ 1.12(5) 1-9(2) 

CH 3 OH 12m 145,766.2 16(0,16)-16(-1,16) E 0.03(2) 2.0(8) 2.4±0.3xl0 14 1.5±0.3xl0 29 ~1.5xl0 -2 

148.111.9 15(0,15)-15(-1,15) E 0.04(2) 3.0(8) 
150,141.6 14(0,13)-14(-1,14) E 0.04(1) 1.9(8) 
151,860.2 13(0,13)-13(-1,13) E 0.05(1) 2.5(8) 

157.048.6 6(0,6)-6(-l,6) E 0.23(1) 2.4(8) 

OH3OH AVERAGE: 5.7±0.7xl0 14 1.3±0.3xl0 29 ~1.3xl0 -2 

CH3CN 12m 147,129.2 8(4)-7(4) 0.019(1) 1.6(8) 2.6±0.3xl0 12 1.7±0.3xl0 27 ~1.7xl0 -4 

147,149.0 8(3)-7(3) 0.054(1) 2.8(8) 

147.163.2 8(2)-7(2) 0.034(1) 3.2(8) 

147.171.7 8(1)-7(1) 0.049(1) 2.4(8) 
147,174.5 8(0)-7(0) 0.051(1) 2.4(8) 

(CH 2 OH) 2 12m 147,131.8 15(7,*,0)-14(7,*,0) <0.01 

Comet T7/LINEAR 

CH 3 OH 12rn 157,246.1 6(0,6)-6(-l,6) E 0.087(1) 1.9(10) 2.2±0.3xl0 Ia 2.0±0.3xl0 2Y ~1.5xlO~ 2c 

157,270.9 1(0,1)-1(-1,1) E 0.051(1) 1,9(10) 

157.272.3 3(0,3)-3(-l,3) E 0.101(1) 1,9(10) 
157,276.0 2(0,2)-2(-l,2) E 0.058(1) 1,9(10) 

a For the BIMA data, In ne represents AJ and are in units of Jy beam - . For the 12m observations, In ne represents T^ in units of K. 
b Comct Hale-Bopp: Q(_ff 2 O) = 10 31 s -1 (Crovisier ct al. 2004) 
c Comct T7/LINEAR: Q(H 2 0)=1.3X 10 29 s _1 (Milam ct al. 2006) 
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Fig. 1. — CH3OH lines detected toward Comet Hale-Bopp with the BIMA array, (a) Emis- 
sion contours from the 2(0,2)-l(0,l) E transition of CH 3 OH at 96.741 GHz. Contours indicate 
the location of the 2(0,2)-l(0,l) E CH 3 OH emission. The contour levels are -0.2, 0.4, 0.6, 
0.8, 1.0, and 1.2 Jy beam -1 . The synthesized beam of 10'/3 x 8'/5 is indicated at the bottom 
left corner, (b) CH3OH spectrum toward Comet Hale-Bopp (Hanning smoothed over three 
channels). The rms noise level is ~0.1 Jy beam -1 (indicated by the vertical bar at the left 
of the spectrum). The spectral line labels correspond to the rest frequency located at the 
top left of the spectral window for a velocity of km s -1 . The dashed line is centered on 
this velocity. The blue trace shows a model fit (see text), (c) Emission contours from the 
3(l,3)-4(0,4) A + transition of CH 3 OH at 107.014 GHz. Contours indicate the location of the 
CH3OH emission averaged between -2 km s -1 and 2 km s -1 . The contour levels are -0.05, 
0.20, 0.30, 0.40, 0.50, and 0.60 Jy beam -1 . The synthesized beam of 10"8 x 7"2 is indicated 
at the bottom left corner, (d) CH 3 OH spectrum toward Comet Hale-Bopp. The rms noise 
level is ~0.1 Jy beam -1 (indicated by the vertical bar at the left of the spectrum). The 
spectral line label is the same as in (b) and the blue trace shows a model fit (see text). 

Fig. 2. — CH3OH lines detected toward Comet Hale-Bopp with the 12m telescope. The 
backends employed for the presented observations were either filter banks with 500 kHz 
resolution or a millimeter autocorrelator (MAC) with a resolution of 781 kHz. The spectral 
line labels correspond to the rest frequency located at the top left of the spectral window for 
a velocity of km s -1 (dashed line). The blue traces show a model fit (see text). 

Fig. 3. — CH3OH lines detected toward Comet LINEAR with the 12m telescope. Spectral 
resolution is 500 kHz and the spectral line labels are similar to Figure 2. The blue traces 
show a model fit (see text). 

Fig. 4. — CH 3 CN lines detected toward Comet Hale-Bopp with the 12m telescope. Spectral 
resolution is 781 kHz and the spectral line labels are similar to Figure 2. The blue trace 
shows the model fit (see text). 

Fig. 5. — CH3OH transitions from Comet LINEAR and the model fit to the data (blue trace) 
from the BIMA array data. 

Fig. 6. — CH3OH Haser model fit to the BIMA array data, (a) The emission contours from 
the 3(l,3)-4(0,4) A+ transition of CH 3 OH at 107.014 GHz. (b) The Haser model predictions 
to these data assuming a scale length of r~10" and then convolved with the synthesized 
beam (shown at the bottom left), (c) The residuals obtained by subtracting the model from 
the data. 
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